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EXECUTIVE  SUMMARY 


This  report  describes  modifications  made  to  the  USARIEM  model  to 
extend  its  applicability  to  high  terrestrial  altitude  environments.  Current 
deployments  of  U.S.  forces  to  high  altitude  regions  in  and  around  Afghanistan 
provided  the  immediate  impetus  for  this  effort.  Primary  focus  has  been  placed  on 
elevations  ranging  from  0  (sea  level)  and  4000  meters  (13,123  ft).  The  approach 
draws  exclusively  on  previously  published  work,  and  is  necessarily  based  on  a 
combination  of  rational  and  empirically  derived  relationships  compatible  with  the 
USARIEM  model's  computational  structure.  The  modifying  algorithms  use 
atmospheric  pressure  and  are  applied  to  the  convective  and  evaporative  heat 
transfer  components  of  the  USARIEM  model,  specifically  the  thermal  resistance 
[It  ]  and  maximum  evaporative  power  of  the  environment  [Emax  ]  algorithm 
derivations.  These  two  modifications  are  minimally  invasive  to  the  USARIEM 
model  computational  engine  and  should  be  applicable  to  hyperbaric  as  well  as 
more  drastic  hypobaric  environments. 

The  additional  input  parameter  required  for  Heat  Stress  Monitor  (HSM) 
instantiations  of  the  USARIEM  model  is  the  barometric  pressure  measured  by  a 
new  on-board  sensor  calibrated  in  atmospheres  (1  ATA=  760  Torr).  For  non-HSM 
instantiations,  barometric  pressure  is  computed  from  terrain  elevation; 

Patm  =(  1  -  2.5577  •  10'^  •  Z  f  where  Patm  is  pressure  in  atmospheres  and 
Z  is  terrain  elevation  in  meters. 

Representation  of  altitude  effects  on  convective  heat  transfer  is  accomplished 
by  modifying  the  computation  of  total  do,  It  (Gonzalez,  1988),  to  include  an 
atmospheric  pressure  term  that  corrects  for  changes  in  air  density.  Measured  or 
computed  Patm,  in  atmospheres,  is  inserted  as  a  simple  multiplier  for  VEtf  in  the 
expression  to  compute  the  altitude/pressure-sensitive  total  do;  ItA  =  Itc  *  (Patm  * 
VEff  )'‘''‘^  where  VEff  is  the  effective  air  velocity  in  m*s"^  and  Itc  and  Itvc  are  the  do 
and  do  velocity  coefficients  respectively  for  the  selected  uniform 

Representation  of  altitude  effects  on  evaporative  heat  transfer  is  accomplished 
by  modifying  the  computation  of  maximum  evaporative  capacity,  Emax ,  to  include 
an  atmospheric  pressure  term.  Measured  or  computed  Patm,  in  atmospheres,  is 
used  to  produce  a  correction  factor  for  the  Lewis  Relation  (LR=  2.2  °C/Torr) 
embedded  in  the  original  Emax  expression  [Emax  =LR  *6.45  •  Ad  •  Cevap  *  (Pssk  - 
Pair)]-  The  altitude/pressure-sensitive  Emax  is  computed  as;  EmaxA  =  Patm’°‘^^ ' 
Emax  (Nishi  and  Gagge,  1977). 

The  result  is  that  dry  thermal  resistance  increases  with  altitude,  as  does 
Emax.  Thus  if  all  other  variables  are  constant,  dry  heat  loss  decreases  with 
altitude  but  total  skin  evaporation  would  increase.  Dry  skin  in  cool  environments 
is  a  little  better  insulated  (with  added  clothing)  and  protected  from  cold  injury  with 
increasing  altitude.  In  hot  and  or  sweating  conditions  where  generally  dry  heat 
loss  is  small  in  comparison  to  evaporative  heat  loss,  evaporative  capacity 
increases  diminishing  heat  stress  and  or  extends  the  person’s  safe  working  time 
for  the  environment  and  activity. 
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INTRODUCTION 


OVERVIEW 

The  USARIEM  Heat  Strain  model  has  been  widely  used  to  develop  tailored 
heat  stress  management  guidelines  for  the  prevention  of  heat  injury  in  military 
training  and  operational  settings  (FM  21-10,  deployment  manuals).  Recently, 
the  USARIEM  model  has  been  implemented  in  a  hand-held  Heat  Stress  Monitor 
(HSM)  as  well  as  the  command  and  control  oriented  Operational  Medicine 
Environmental  Grid  Applications  (OMEGA)  test  bed  system  (5,  6).  The  model's 
principal  output  parameters  are  optimal  work/rest  schedules  and  maximum  safe 
work  time  limits  based  on  predicted  body  temperature  changes,  and  hourly 
drinking  water  needed  to  replace  predicted  fluid  losses  from  sweating.  The 
primary  inputs  are  environment  (air  temperature,  humidity,  wind  speed,  and  solar 
load),  clothing  characteristics  (insulation  and  vapor  permeability),  metabolic  rate, 
and  acclimatization  status. 


STATEMENT  OF  THE  PROBLEM 

The  USARIEM  model  is  fundamentally  an  empirical,  operational  model 
applied  to  temperature  to  heat  stress  climatic  zones.  It  is  based  on  the  results  of 
many  years  of  laboratory  and  field  studies  with  human  volunteers  over  a  wide 
range  of  hot  environments,  but  at  locations  generally  within  a  few  hundred 
meters  of  sea  level.  Consequently,  the  wet  and  dry  heat  transfer  functions 
represented  in  the  model  do  not  include  consideration  of  the  effects  of  high 
terrestrial  altitude  on  heat  strain.  The  military  response  to  global  terrorist 
activities  will  likely  continue  to  involve  deployments  to  geographic  regions  that  lie 
significantly  above  sea  level  and,  because  such  deployments  may  present  the 
added  heat  injury  risks  associated  with  protective  clothing  encapsulation,  there  is 
a  need  to  extend  the  USARIEM  model's  applicability  at  altitude. 

OBJECTIVES 

The  primary  objective  of  this  report  is  to  document  modifications  to  the 
USARIEM  heat  strain  model  that  are  intended  to  extend  its  applicability  to  heat 
stress  conditions  in  high  terrestrial  altitude  environments.  A  secondary  objective 
is  to  illustrate  the  results  of  those  modifications  on  the  model’s  output  in  a  context 
relevant  to  conventional  and  chemical  protective  clothing  scenarios.  Scope  is 
limited  to  methods  and  tools,  and  does  not  include  a  rigorous  sensitivity  analysis 
that  will  be  the  focus  of  a  subsequent  report. 

APPROACH 

The  scientific  basis  for  the  modifications  described  here  draws  exclusively 
on  previously  published  work.  The  approach  is  further  constrained  to  a 
combination  of  rational  and  empirically  derived  m  athematical  relationships  that 
are  compatible  with  the  USARIEM  model's  existing  computational  structure.  The 
resulting  modifications  are  focused  on  a  practical  representation  of  ambient 
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atmospheric  pressure  effects  on  the  convective  and  evaporative  heat  transfer 
functions  in  the  model,  specifically  mathematical  algorithms  affecting  the  clothing 
thermal  insulative  properties  (It  )  and  maximum  evaporative  power  of  the 
environment  (Emax  )■  The  required  additional  parameter  for  the  model's 
environmental  input  list  is  barometric  pressure  or  altitude  (hypobaric)  effects 
above  sea  level  altering  the  above  two  factors. 


METHODS 

DEVELOPMENT  PROCEDURES  AND  TOOLS 
Spreadsheet  Version 

A  preliminary  evaluation  of  the  behavior  of  the  modifications  to  the  It  and 
Emax  functions  was  accomplished  using  a  spreadsheet  instantiation  of  the 
USARIEM  model  (Matthew  et  al.,  unpublished).  This  spreadsheet  model  version 
provides  minute-to-minute  time  series  outputs  of  intermediate  calculated 
parameters  and  body  core  temperature  as  well  as  values  for  the  modified  It  and 
Emax  parameters  during  pre-defined  work/rest  activity  scenarios.  This  allows  a 
graphical  view  of  the  impact  of  altitude  predicted  core  temperature  responses 
across  the  relevant  time  domain. 

C  Language  Version 

Source  Code:  The  C  language  instantiation  of  the  USARIEM  Heat  Strain  Model 
is  based  on  transformed  source  code  origianlly  derived  primarily  from  the 
MERCURY/OMEGA  project’s  Heat  Strain  Decision  Aid  (HSDA)  module 
(Matthew,  2000).  Three  source  code  files  are  compiled  and  linked  during  the 
build  process:  1)  erf.c  that  computes  the  error  function,  2)  heatcas.c  that 
estimates  heat  casualty  risk  using  final  core  temperature  and  the  computed  error 
function,  and  3)  hsdac.c  that  contains  the  algorithms  comprising  the  main 
computational  engine. 

Revisions:  Recent  modifications  to  the  main  computational  engine  corrected 
original  sweating  rate  computations  and  systematized  handling  of  solar  radiation 
and  estimating  of  the  Wet  Bulb-Globe  Temperature  (WBGT)  index  components 
(Matthew  et  al.,  2001).  The  current  modifications  to  It  and  Emax  which  enable 
model  sensitivity  to  altitude  are  implemented  in  the  most  recent  version  of 
hsdac.c,  hsdacS.c. 

Running  ariem_a:  The  single  executable  file  resulting  from  a  build  of  hsdacS.c, 
erf.c,  and  heatcas.c  is  called  ariem_a.exe.  The  executable  reads  an  input  file, 
and  automatically  saves  the  results  in  an  output  file. 

Input:  The  input  file,  ariem_a.in,  shown  in  Figure  1,  is  an  eighteen  line  ASCII 
text  file  with  the  numeric  input  values  on  the  left  and  the  identifying  variable,  units 
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of  measure,  and  comments  on  the  left.  The  user  edits  the  left  hand  column  with 
a  text  editor  such  as  Notepad  ®,  then  saves  that  file  for  each  run. 


Figure  1.  Example  of  the  ariem  a.in  file. 


25,56 

!  db  (C) 

90,0 

!  rh  (%) 

1.00 

!  wind  speed  (m/s) 

25,56 

!  mean  radiant  temp  (C) 

4000.0 

!  Terrestrial  altitude  (m) 

425,00 

!  work  rate  ( W)  [  1 05,  1 50,  250,  425,  600] 

0.00 

!  work  external  (W) 

12.00 

!  acclimatization  (days)  [0  -  99] 

1.24 

!  dehydration  factor  [0,  1.24,  2.5,  4.0,  6.0] 

39.0 

!  max  work  temp  (C)  [39.0  historically] 

38,50 

!  max  work  rest  cycle  temp  (C)  [38.5  historically] 

36.50 

!  skin  temp  (C)  [fixed] 

1.08 

!  clothing  parameters:  do  hw_bdu=1.08 

0.47 

!  iclo  hw_bdu=0.47 

-0.27 

!  gc  hw_bdu=  -0.27 

0.41 

!  gi  hw_bdu=  0.4 1 

176.0 

!  soldier  height  (cm)  [152-218] 

70.0 

!  soldier  wieght  (kg)  [50  -  120] 

Output  Parameter  Definitions:  The  output  parameters  are  as  follows; 

MaxWk[min]  =  Maximum  safe  one-time  continuous  work  exposure 
(minutes)-  assumes  extended  recovery  time  in  cool  environment 

MxWtr  [qt/hr]  =  Hourly  drinking  water  required  to  balance  sweat  losses 
during  one-time  max  work  exposure 

W/R  [min]  =  Minutes  of  work  per  hour  allowed  (assuming  balance  of  the 
hour  is  spent  at  rest) 

0  output  means  no  sustainable  work  rest  cycle  is  possible 
60  output  means  work/rest  cycles  not  needed 

CyWtr  [qt/hr]  =  Hourly  drinking  water  required  to  balance  sweat  losses 
during  work/rest  cycle 

0  output  means  work/rest  cycle  not  applicable  -(either  0  or  60 
minutes  output  in  W/R  [min]  ) 
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RstWtr  [qt/hr]  =  Hourly  drinking  water  required  to  balance  sweat  losses 
while  at  rest 

RstWtrSh  [qt/hr]  =  Hourly  drinking  water  required  to  balance  sweat 
losses  while  at  rest  in  the  shade 

Patm  =  Computed  local  pressure  in  atmospheres  (1ATA  =  760Torr) 

Heat  Cas  [%]  =  Heat  casualty  expectation  if  no  action  taken  to  implement 
work/rest  cycle  or  max  work  time  limit 

Tdew  [°C]  =  Dew  point  temperature 

Twba[°C]  =  Aspirated  or  psychrometric  wet  bulb  temperature 

Twbn[°C]  =  Estimated  natural  wet  bulb  as  used  in  computation  of  the  Wet 
Bulb  Globe  Temperature  (WBGT)  index 

Tgrc]  =  Estimated  Black  Globe  temperature  as  used  in  computation  of 
the  Wet  Bulb  Globe  Temperature  (WBGT)  index 

WBGT[°C]  =  0.1  *  Ta  (from  input)  +  0.2  *  Tg[°C]  +  0.7*Twbn[°C] 

MWkRec[min]  =  Estimated  recovery  time  in  minutes  (in  work 
environment)  from  maximum  safe  one-time  continuous  work  exposure,  as 
follows: 

(1 )  0  minutes  output  means  no  cooling  possible:  Recovery  not 
possible  in  that  environment 

(2)  >300  minutes  means  cooling  rate  too  slow:  Recovery  not 
possible  in  that  environment 

MWkRecS[min]  =  Estimated  recovery  time  in  minutes  (in  work 
environment  but  in  the  shade)  from  maximum  safe  one-time  continuous 
work  exposure,  as  follows: 

(1)  0  output  means  no  cooling  possible:  Recovery  not  possible 
in  that  environment 

(2)  >300  minutes  means  cooling  rate  too  slow:  Recovery  not 
possible  in  that  environment 

Output  File;  The  output  file,  ariem_a.out,  is  a  two  line,  tab  delimited  ASCII  text 
file  shown  in  Figure  2.  The  first  line  contains  the  parameter  identifier  text  and  the 
second  line  contains  the  output  values. 
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Figure  2.  Example  of  the  ariem  a.out  file. 


MaxWklmin]  MxWtr  [qt/hr] 

132  0.9 

W/R  [min] 

32 

CyWlr  [qt/hr] 

0.7 

RstWlr  [qVlir] 
0.4 

RslWtrSh  [ql/hr] 

0.2 

Heat  Cas  [%]  . 

24.4  . 

.  'I’dewfC]  Twba[CJ 

.  23.8  24.2 

Twbn[C] 

25.1 

■rg[Cl  WBGT[C] 
46.6  29.5 

MWkRcc[min] 

33 

MWkRccS[minJ  Patm 

22  0.57 

Because  the  output  file  is  overwritten  each  time  the  program  is  run,  it  is 
necessary  to  save  contents  between  runs.  The  horizontal  format  is  intended  to 
facilitate  a  “cut  and  paste”  operation  into  a  standard  spreadsheet  file  to 
accumulate  multiple  model  runs  for  analysis. 


FUNCTION  MODIFICATIONS 

Terrestrial  Altitude  and  Barometric  Pressure 

Altitude-related  modifications  to  the  model  are  based  on  barometric 
pressure.  Our  modifications  assume  local  barometric  pressure  is  not  known  and 
is  approximated  from  known  terrain  elevation  (1); 

Patm  =(1  -  2.5577  •  10'^  •  Z  )^2559  Atmospheres  (Eq.1) 

Where, 

Patm  =  pressure,  atmospheres 
Z  =  terrain  elevation,  meters 

This  addition  to  the  source  code  is  implemented  as  a  defined  function,  PATM,  in 
the  header  section  of  hsdacS.c.  Terrestrial  altitude  in  meters  is  read  from  the 
input  file  and  Patm  is  computed  as  described  above.  It  should  be  noted  that  this 
function  is  unnecessary  in  HSM  and  OMEGA  applications  because  automated 
measurements  of  barometric  pressure  would  allow  direct  input  of  Patm  rather  than 
Z. 


Convective  Heat  Transfer 

Representation  of  altitude  effects  on  convective  heat  transfer  is 
accomplished  by  modifying  the  computation  of  total  do,  k,  to  include  an 
atmospheric  pressure  term  that  corrects  for  changes  in  air  density  (1,3,8). 
Measured  or  computed  Patm,  in  atmospheres,  is  inserted  as  a  simple  multiplier  for 
VEff  in  the  expression  to  compute  the  altitude/pressure-sensitive  total  do: 
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ItA  -  Itc  *  (Patm  •  Veff  ) 


do 


(Eq.  2) 


Itvc 


Where 

P  atm  =  pressure,  atmospheres 

ItA  =  air  velocity  and  pressure  adjusted  do  for  the  selected  uniform 
Veff  =  effective  air  velocity,  m*s'^ 

Itc  =  do  at  1  m*s"^  for  the  selected  uniform 

Itvc  =  air  velocity  modifier  for  do  for  the  selected  uniform 

Evaporative  Heat  Transfer 

Representation  of  altitude  effects  on  evaporative  heat  transfer  is 
accomplished  by  modifying  the  computation  of  maximum  evaporative  capacity, 
Emax-  Specifically,  an  atmospheric  pressure  term  is  added  to  adjust  the  value  of 
the  Lewis  Relation  (LR)  embedded  in  the  original  E^ax  expression  (3,8). 
Measured  or  computed  Patm,  in  atmospheres,  is  used  to  produce  a  correction 
factor  for  LR  as  previously  described  (1),  and  the  resulting  altitude/pressure- 
sensitive  Emax  is  computed  as: 

EmaxA  “  P atm  -0-45  .  Lp  .  0  45  ,  Watts  ( Eq.  3) 

where, 

EmaxA  =  atmospheric  pressure  sensitive  Emax,  Watts 
Patm  =  pressure,  atmospheres 
LR  =  Lewis  Relation,  2.2  °C  •  Torr’’' 

Ad  =  Dubois  body  surface  area,  m^ 

Cevap  =  water  vapor  permeability 

Pssk  =  saturation  vapor  pressure  at  Tskin ,  Torr 

Pair  =  ambient  air  vapor  pressure,  Torr 

RESULTS 

THE  SPREADSHEET  IMPLEMENTATION 
Body  Core  Temperature  Profiles 


Figure  3  shows  spreadsheet-predicted  body  core  temperature  profiles  for  heat 
stress  exposure  at  various  terrestrial  elevations  wearing  the  battledress  uniform 
(BDU). 
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Figure  3.  Spreadsheet  predicted  Tr©  during  Work/Rest  cycles  in  BDU  at  425 
Watts. 

Time  series  input  scenario  begins  with  a  ten  minute  rest  followed  by  five 
consecutive  cycles  of  30  minutes  of  work  (M=  425  Watts  )  and  30  minutes  of  rest 
(M=  105  Watts).  Air  temperature  is  40  °C,  RH  is  20%,  wind  speed  is  2  m/s,  and 
mean  radiant  temperature  is  60  °C 


Figure  4  shows  spreadsheet-predicted  body  core  temperature  profiles  for  heat 
stress  exposure  at  various  terrestrial  elevations  while  wearing  the  new  joint 
services  light  integrated  suit  technology  (JSLIST)  uniform  over  the  hot  weather 
BDU.  Time  series  input  scenario  again  begins  with  a  ten  minute  rest  followed  by 
five  consecutive  cycles  of  30  minutes  of  work  (M=  405  Watts  )  and  30  minutes  of 
rest  (M=  105  Watts).  Air  temperature  in  this  scene  is  30  °C,  RH  is  40%,  wind 
speed  is  2  m/s,  and  mean  radiant  temperature  is  60  °C. 
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Figure  4.  Spreadsheet  predicted  Tre  during  Work/Rest  cycles  in  JSLIST 
over  BDU  at  425  Watts. 


Modulated  by  time  delay  algorithms  which  change  whenever  there  is  a  transition 
to  a  different  value  for  metabolic  rate,  the  core  temperature  profile  is 
fundamentally  driven  by  the  model’s  internally  predicted  final  equilibrium  rectal 
temperature  parameter,  Tre.f.  It  represents  the  theoretical  final  steady  state  rectal 
temperature  for  a  static  condition  set.  For  the  BDU  scenario  the  sea  level  and 
4000  meter  Tre, f  values  were  37.55  and  37.46  °C  respectively  at  105  Watts  and 
39.07  and  38.82  C  at  425  Watts.  For  the  JSLIST  scenario  the  sea  level  and 
4000  meter  Tre,f  values  were  37.75  and  37.64  °C  respectively  at  105  Watts  and 
40.40  and  39.92  °C  at  425  Watts.  The  model  output  clearly  suggests  that  the  net 
effect  of  the  high  terrestrial  altitude  is  to  improve  heat  transfer  and  mitigate  heat 
strain  in  both  conventional  (BDU)  and  protective  encapsulation  (JSLIST) 
scenarios. 
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Water  Requirements 

As  in  the  C  language  version  of  the  ARIEM  operational  model,  prediction  of  water 
requirements  in  the  spreadsheet  version  has  no  time  dependency.  It  is  based  on 
prediction  of  an  average  equilibrium  sweat  rate  (Shapiro,  1982)  that  remains 
constant  for  a  particular  environment,  clothing  type,  and  work  load.  Response 
lag  time  algorithms  have  not  been  implemented  and,  consequently,  transitions 
are  abrupt  when  viewed  in  time  series  format.  Figure  5  shows  the  spreadsheet- 
predicted  hourly  drinking  water  requirements  for  the  working  (M=  425  Watts)  and 
resting  (M=105  Watts)  conditions  in  BDU.  Air  temperature  is  40  °C,  RH  is  20%, 
wind  speed  is  2  m  s'\  and  mean  radiant  temperature  is  60  °C. 

Figure  5.  Spreadsheet-predicted  water  requirements  in  BDU  at  425  Watts 
and  at  rest. 
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Figure  6  shows  the  spreadsheet-  predicted  hourly  drinking  water  requirements 
for  the  working  (M=  425  Watts)  and  resting  {M=105  Watts)  conditions  in  JSLIST 
worn  over  BDU.  Air  temperature  is  30  °C,  RH  is  40%,  wind  speed  is  2  m/s,  and 
mean  radiant  temperature  is  60  °C. 
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Figure  6.  Spreadsheet-predicted  water  requirements  in  JSLIST  over  BDU  at 
425  Watts  and  at  rest. 
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Differences  between  the  sea  level  and  4000  meter  water  requirements  are  quite 
small  for  both  the  conventional  (BDU)  and  protective  encapsulation  (JSLIST) 
scenarios. 


THE  C  LANGUAGE  IMPLEMENTATION 
Source  Code 

Source  code  listings  for  this  updated  C  language  version  of  the  USARIEM 
heat  strain  model  are  provided  at  Appendix  A.  The  source  code  listing  for 
hsdacS.c  contains  several  blocks  of  ‘dead  code’  that  have  been  retained  to  help 
document  revision  history.  While  the  core  functionality  of  the  code  is  adequate 
for  the  simplest  analytical  tasks,  it  is  clear  that  a  rigorous,  disciplined  re-write  of 
the  computational  engine  is  in  order.  The  addition  of  a  graphical  user  interface 


11 


(GUI)  and  input/output  facilities  to  automate  multiple  ‘case’  runs  would  enhance 
the  overall  usefulness  of  the  model.  Future  modifications  are  planned  to 
implement  this  aspect  of  the  model. 

Output  Products 

The  executable  C  language  version  of  the  USARIEM  heat  strain  model 
(ariem_a,  compile  date  11/30/2001 )  provides  tabular  output  for  a  single  input 
situation  or  ‘case’.  Output  for  each  ‘case’  includes  maximum  safe  work  time, 
optimal  work/rest  cycle  limit,  associated  water  requirements  as  well  as  a  heat 
casualty  probability  risk  assessment,  and  a  computed  estimate  of  the  ambient 
wet  bulb-globe  temperature  (WBGT)  index.  Table  1  depicts  output  results  for  a 
clothing  configuration  consisting  of  the  JSLIST  uniform  worn  over  the  BDU.  Each 
line  in  the  table  used  the  ariem_a  executable  \A/ith  a  different  input  file  and 
contents  of  the  resulting  output  file  were  ‘pasted’  into  a  summary  spreadsheet 
file.  That  process  was  repeated  for  each  of  the  five  terrestrial  altitudes  within 
each  of  the  three  assumed  metabolic  rates 


Table  1.  The  ariem  a  model  output  sensitivity  to  high  terrestrial  altitude  in 
chemical  protective  posture  at  three  different  work  rates. 

JSLIST  over  BDU  at  Ta=  30  "C,  RH=40%,  Va=2  m/s,  Tmr=60  “C 


M=250  Watts 

Altitude  (m)  MxWk[min] 

CyWtr  [qt/hr] 

Heat  Cas  [%] 

0 

300 

0.8 

0.8 

7.7 

1000 

300 

0.8 

6.4 

2000 

300 

■  0.8 

5.2 

3000 

300 

'  0.8 

4.2 

4000 

300 

0.7 

3.3 

M=350  Watts 

Altitude  (m)  MxWk[min] 

MxWtr  [qt/hr]  W/R  [min] 

CyWtr  [qt/hr] 

Heat  Cas  [%] 

0 

81 

1.1 

24 

0.7 

59.0 

1000 

85 

1  1 

■■  25 

0.7 

52.8 

2000 

91 

1.1  ' 

27 

0.7 

46.4 

3000 

98 

1.0 

28 

0.7 

39.9 

4000 

108 

1.0 

31 

0.7 

33.5 

M=425  Watts 

Altitude  (m)  MxWk[min] 

MxWtr  [qt/hr] 

W/R  [min] 

CyWtr  [qt/hr] 

Heat  Cas  [%] 

0 

58 

1.3 

19 

0.7 

96.1 

1000 

60 

1.3 

19 

0.7 

93.9 

2000 

62 

1.3 

20 

0.7 

90.6 

3000 

64 

1.2 

21 

0,7 

86.1 

4000 

68 

1.2 

23 

0.7 

80.2 
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DISCUSSION 


VERIFICATION  AND  VALIDATION 

Verification  issues  for  the  USARIEM  heat  strain  model  are  somewhat 
complicated  by  its  empirically-derived,  population-based,  development  history.  In 
order  to  provide  a  margin  of  safety  for  those  individuals  whose  response  would 
be  greater  than  the  statistical  average  response,  the  model  was  intentionally 
‘tuned’  to  over-predict  the  final  equilibrium  Tref  by  as  much  as  two  standard 
deviations.  The  resulting  optimal  work/rest  cycles  and  maximum  safe  work  time 
guidance  are  therefore  very  conservative.  While  this  is  a  desirable  and  probably 
necessary  characteristic  for  a  tactical  decision  aid  application,  it  results  in  a 
somewhat  disconcerting  level  of  error  when  averaged  field  or  laboratory 
response  measurements  are  compared  with  model  output.  Previous  work  has 
indicated  that  an  adjustment  to  the  denominator  of  the  Kwork  expression  from  120 
to  225  is  sufficient  to  bring  the  time  dependent  outputs  into  line  with  averaged 
response  measurements  (Gonzalez,  1997;  Cadarette,  1999).  The  ability  to 
switch  between  the  conservative  and  statistical  average  response  formulations 
for  Kwork  depending  on  application  context  (i.e.  either  operational  or 
research/validation)  would  extend  the  global  utility  of  the  model.  Nevertheless,  a 
systematic  revision  of  the  model  to  accurately  predict  mean  Tre  response  coupled 
with  a  statistically  valid  stochastic  representation  of  core  temperature 
probabilities  in  the  high,  or  “risk  tail”  portion  of  the  response  distribution  would 
probably  result  in  a  more  robust  risk  assessment  paradigm.  For  the  sake  of 
consistency  with  a  large  number  of  recent  modeling  analysis  tasks,  the  current 
instantiation  employs  the  original,  conservative  120  value  in  the  Kwork  expression. 

USER  APPLICATIONS 

The  range  of  thermal  strain  applications  for  the  USARIEM  model  is 
considerably  expanded  with  the  addition  of  the  algorithms  to  instantiate 
sensitivity  to  high  terrestrial  altitude  following  basic  heat  transfer  principles 
outlined  by  previous  work  (Nishi  and  Gagge,  1977;  Gonzalez,  1988).  While 
these  algorithms  do  not  consider  non-thermal  health  risks  associated  with  living 
or  working  at  high  terrestrial  altitudes,  they  do  extend  the  heat  strain  model 
working  domain  to  terrestrial  elevations  that  are  consistent  with  emerging  geo¬ 
political  deployment  scenarios.  The  additional  improvements  and  additions  to 
the  current  ARIEM  heat  stress  model  is  ideally  applicable  to  various  Objective 
Force  Warrior  requirements  and  Land  Warrior  clothing  systems  planned  for  the 
future.  The  implementation  of  the  model  in  C  language  options  is  also  a  marked 
improvement  from  the  previous  versions. 
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CONCLUSIONS 


The  simple  ASCII  text  input/output  file  interface  to  the  ariem_a  model  severely 
limits  its  utility  for  anything  more  than  the  simplest  analytical  tasks.  There  is  a 
need  to  design  and  implement  a  Graphical  User  Interface  (GUI)  and  extensions 
that  allow  a  batch  mode  capability  that  will  automate  preparation  of  a  single, 
larger  input  file  consisting  of  user-specified  multiple  domain  ranges  and  step 
intervals 
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APPENDIX  A 


heatcas.c 

double  erf (double  XX) ; 


float  heatcas (float  Tmxcore) 

/*********************************************************************** 
■k  -k  -k  -k  -k 

computes  heat  strain  casualty 

■kkkk-k-k-k-k-k-kkk-kk-k-kkk-kk-kk-k-k-k-k-k-k-k-k-k-k'k-kkkkkk-kk'k-k-kk-kkick-kkk-kk-k-kkk'k-kick-kk-kkk-k-kkkk 
■k  k  k  k  ^ 

{ 

static  float  Tmean  =  39. 5F,  sigma  =  0.5198F,  sqrt2  =  1.4142135F; 
float  heat_cas,  arg; 

arg  =  (Tmxcore  -  Tmean)  /  (sqrt2  *  sigma) ; 
heat_cas  =  (float) erf (arg)  /  2.0F  +  0.5F; 

heat_cas  *=  100. OF;  /*  =>  %  */ 

return (  heat_cas  ) ; 

} 


Erf.c 


#include  <math.h> 


double  erf (double  XX) 

/********************************************************************** 
C  Computes  the  value  of  the  error  function  at  the  point  XX 
Ref:  C.  Hastings,  Approximations  for  Digital  Computers, 

Princeton  Univ.  Press,  Princeton,  NJ,  1955 
Referenced  by:  DSDIST,  XYPVPR 
External  References :  None 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk^ 

{ 


double  ERFCTN,  SIGN,  X,  X2 ,  X3 ,  X4 ,  X5 ,  X6 ,  Y; 

static  double  A1=0 . 0705230784 ,  A2=0 . 0422 820123 ,  A3=0 . 0092705272 , 
A4  =  0 . 0001520143 ,  A5  =  0 . 0002765672 ,  A6  =  0 . 0 0 0 04 3 063 8  , 
PI=3 . 14159265,  XMIN=l.E-5,  XMAX=5 . ; 


X=XX  ; 

if  (X  ==  0  .  ) 
return (0.0) ; 


SIGN  =  X  /  fabs(X); 
X  =  fabs (X) ; 


if  (X  <  XMIN) 

ERFCTN  =2.  *  X  /  sqrt(PI); 

else  if  (X  >  XMAX) 
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ERFCTN=1.0; 
else  { 

X2=X*X; 

X3=X2*X; 

X4=X3*X; 

X5=X4*X; 

X6=X5*X; 

Y=1 . +A1*X; 

Y=Y+A2*X2 ; 

Y=Y+A3*X3 ; 

Y=Y+A4*X4 ; 

Y=Y+A5*X5; 

Y=Y+A6*X6; 

Y  =  pow ( Y ,  16.0); 

Y  =  1.  /  Y; 

ERFCTN  =  1 .  -  Y ; 

} 

X  =  SIGN  *  X; 

ERFCTN  =  SIGN  *  ERFCTN; 

return (  ERFCTN  ) ; 

} 

hsdacS.c 

/********************************************************************** 
program:  ariem_a;  This  file  is  HSDACS.C  : Computational  Engine  Heat 
Strain  Decision  Aid  C,  version:  8 
created:  October  2001 

History , USARIEM  model  source  code: 

BASIC  Language  engine-  L.A.  Stroschein,  USARIEM: 

C  Language  PC  version  of  firmware  implementation  for  Heat  Stress  Monitor 
(HSM)  -  K.  Honeyager,  Under  Contract,  SwRI 

C  Language  functional  extensions  for  MERCURY/OMEGA  -  J.R.  Furlong,  Under 
Contract,  SAIC 

C  Language  engine  modifications  for  high  terrestrial  altitude  -  W. 
Matthew  &  L.  Berglund,  USARIEM 

REFERENCES 

1)  Pandolf,  K.B.,  L.A.  Strochein,  L.L.  Drolet,  R.R.  Gonzalez,  M.N. 

Sawka,  "Predictive  Modeling  of  Physiological  Responses  and  Human 
Performance  in  the  Heat,"  Comp.  Biol.  &  Med.,  6:319-329,  1986. 

#include  <stdio.h> 

#include  <string.h> 

#include  <math.h> 


ttdefine  data 
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#define  code 
#define  xdata 
#define  idata 


#define  fin  "ARIEM_a.in" 

#define  font  " ARIEM_a . out " 

#define  LBUF  128 

#define  DUBOIS (ht , wt )  0.007184F  *  (f loat ) (pow (ht,  0.725) 
0.425)) 

ttdefine  PATM(alt_m)  (float) pow ( (1 . 0  -0.000025577  *  alt_m 

#define  LEWIS_C0R (patm)  (float) (pow(patm,  -0.45)) 

ttdefine  SUCCESS  1 
#define  FAILURE  0 
ftdefine  STEFAN_BOLTZMANN  5.67e-8 

float  heatcas ( float  Tmxcore) ; 

float  tmtry; 

float  Ereqs ; 

float  Emaxs; 


typedef  struct  inputs  { 
float  Ta; 
float  rh; 
float  ws; 
float  mrt; 

float  alt_m; 
float  workrate; 
float  workextn; 
float  acclim; 
float  dehyd; 

float  Tmxwork,  Tmxcyc,  Tnixmetrate,  Tskin; 
float  do,  iclo,  gc,  gi; 
float  ht,  wt ; 

}  INPUTS ; 

INPUTS  hsm_in; 

FILE  *stream; 

char  buf f er [LBUF+1] ; 


typedef  struct  algo__variables  { 


float 

TrefWK; 

float 

TrefRY; 

float 

TrefAWK; 

float 

TrefARY; 

float 

WtrWK; 

float 

WtrRY; 

float 

WtrRYSh; 

float 

CP; 

*  pow(wt, 
, 5 .2559) 
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float  CPSh; 
float  TDWK; 
float  TreoCWK; 
float  KWK; 
float  DTreWK; 
float  TDRY; 
float  KRY; 
float  DTreRY; 
float  TmCWK; 
float  TmCRY; 

float  TreoCRY; 
float  Tskin; 

}  HSMALGO_STRUCT ; 

/* 

*  Function  prototypes  for  Heat  Stress  Monitor  functions . 

*/ 

void  CalcBaseValues  {HSMALGO_STRUCT  *) ; 

void  DoTheRest  (HSMALGO_STRUCT  *) ; 

void  ComputeAspiratedWetBulb  (float  *)  ; 

void  GetSkinTemp  (HSMALGO_STRUCT  *) ; 

void  ComputeDewPoint ()  ; 

int  ComputeNaturalWetBulb ( ) ; 

int  HSPrediction (void) ; 

int  ComputeGlobe  (void) ; 

float  SVP  (float  T) ; 

float  Veff  (float) ; 

float  It (float , float) ; 

float  Cevap  (float) ; 

float  U  (float) ; 

float  Hrc  (float,  HSMALGO_STRUCT  *,  float  Asold) ; 
float  Ereq  (float,  float,  float,  float); 

float  Emax  (float  patm,  float,  float,  HSMALGO_STRUCT  *,  float  Asold); 

float  Tref  (float,  float,  float,  float,  float,  float); 

float  DTref  (float,  float) ; 

float  WTR  (float,  float,  float  Asold); 

float  MxWK  (float,  float,  float,  float,  float); 

float  MxWRKRY  (float,  float,  f loat , float) ; 

float  TreRY  (float  t,  float  TreoCRY,  float  DTreRY,  float  KRY,  float 
TDRY) ; 

float  TreWK  (float  t,  float  TreoCWK,  float  DTreWK,  float  KWK,  float 
TDWK) ; 

float  round  (float,  unsigned  char) ; 

float  f_wetbulb (float  hconv,  float  hevap,  float  Ta,  float  Trk4 ,  float 
Pdp, 

float  Twb) ; 


/* 

*  Structure  definitions. 

*/ 

typedef  struct  clothing_struct  { 
char  *text; 
float  itc; 
float  itvc; 
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float  imc; 
float  imvc; 

}  clothing_table ; 

typedef  struct  clothing_types  { 
char  *text; 

clothing_table  *clo_lst; 

}  clothing; 

/*typedef  struct  work_struct  { 
char  *text; 
float  metabolic_rate ; 

}  work_table ; */ 

struct  select  { 

unsigned  char  clothing_type_idx; 
unsigned  char  clothing_idx; 
/*unsigned  char  work_rate_idx; */ 

in¬ 


struct  hsm_data  { 
float  globe; 
float  wbgt ; 
float  wet_bulb; 
float  awet_bulb; 

float  dry_bulb; 
float  wind_speed; 
float  relhum; 
float  mrt; 

float  alt_m; 
float  dewpoint; 

float  Lcor; 

}; 


/* 

*  Define  the  MISCELLANEOUS  clothing  option  list. 
*/ 

clothing_table  code  misc__list  [1]  =  ( 

"MERCURY  INPUT  ",  -99.0,  -99.0,  -99.0, 

in¬ 


clothing  code  clothing_list [1]  =  { 

"MISCELLANEOUS  ",  misc_list, 

in- 


/ *work__table  code 
"LEVEL  1 
"LEVEL  2 
"LEVEL  3 
"LEVEL  4 
"LEVEL  5 
NULL, 

/*};*/ 


work_list  []  =  { 

",  105.0, 

",  150.0, 

",  250.0, 
",  425.0, 
",  600.0, 
0 . 0*/ 


-99.0, 
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/ 


*  structure  containing  menu  selections  used  in  the  calculation  of 

*  work/rest  and  water  requirements. 

*/ 

xdata  struct  select  input_selections ; 
struct  hsm_data  measurements; 

/*  global  variables  */ 
xdata  float  MRT; 


xdata  float  patm; 

/*  Heat  stress  prediction  results  */ 

xdata  int  work_rest; 

xdata  float  water_req; 

xdata  int  max_work; 

xdata  float  max_water; 

xdata  float  rest_water; 

xdata  float  rest  water  sh; 


xdata  float  maxwk_rcy; 

xdata  float  maxwk_rcy_sh; 
xdata  float  heat  cas; 


main  (int  argc,  char  **argv) 


{ 

/*int  argi;  */ 

/*unsigned  char  string[80];  */ 
clothing_table  *cl_table; 


/*  read  ariem_a.in  */ 
stream  =  fopen{fin,  "r"); 
if  (stream  ==  NULL)  { 

strcpy (buf f er ,  fin)  ; 
printf("Can  not  open  %s\n",  buffer); 
exit { 1 ) ; 

} 


fgets (buffer,  LBUF,  stream) ; 
sscanf (buffer,  "%f",  &hsm_in.Ta); 
fgets (buf fer,  LBUF,  stream); 
sscanf (buffer,  "%f",  &hsm_in.rh); 
fgets (buf fer ,  LBUF,  stream); 
sscanf  (buffer,  "%f",  &;hsm_in.ws); 
fgets (buf fer ,  LBUF,  stream); 

sscanf (buffer,  "%f",  &hsm_in.mrt) ; 

fgets (buffer,  LBUF,  stream) ; 
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sscanf (buffer, 


%f",  &hsm_in.alt_m); 

fgets (buf f er ,  LBUF,  stream); 

sscanf  (buffer ,  "%f",  &:hsm_in.workrate); 

fgets (buf fer,  LBUF,  stream); 

sscanf  (buffer,  "%f",  &lism_in.workextn); 

fgets (buf fer ,  LBUF,  stream); 

sscanf (buffer,  "%f",  &hsm_in.acclim); 

fgets (buf fer ,  LBUF,  stream); 

sscanf  (buffer,  "%f",  &:]ism_in  .  deliyd)  ; 

fgets (buf fer,  LBUF,  stream); 

sscanf  (buffer,  "%f",  &]ism__in .  Tmxwork)  ; 

fgets (buffer,  LBUF,  stream) ; 

sscanf  (buffer,  "%f",  &]ism__in.Tmxcyc); 

fgets (buf fer ,  LBUF,  stream); 

sscanf  (buffer,  "%f",  &lism_in.Tmxmetrate); 

fgets (buf fer ,  LBUF,  stream); 

sscanf  (buf  fer ,  "%f",  &lism_in.Tskin); 

fgets (buf fer ,  LBUF,  stream); 

sscanf (buffer,  "%f",  &hsm_in.clo); 

fgets (buf fer ,  LBUF,  stream); 

sscanf  (buffer,  "%f",  &lism_in.iclo); 

fgets (buf fer ,  LBUF,  stream); 

sscanf  (buffer,  "%f",  &lism_in.gc); 

fgets (buf fer ,  LBUF,  stream); 

sscanf (buffer,  "%f",  &hsm_in.gi); 

fgets (buf fer ,  LBUF,  stream); 

sscanf  (buffer,  "%f",  &lism_in.lit)  ; 

fgets (buf fer,  LBUF,  stream); 

sscanf  (buf fer ,  "%f",  &lism_in . wt )  ; 

fclose (stream) ; 

measurements . globe  =  42.8896F; 
measurements  .  dry_bulb  =  ]ism_in.Ta; 
measurements . relhum  =  hsm_in.rh; 
measurements . wind_speed  =  hsm_in.ws; 

measurements  . mrt  =  lism_in.mrt;  /*  [C]  */ 

MRT  =  ]ism_in.mrt;  /*  global  MRT  */ 

patm  =  PATM (hsm_in. alt_m) ; 

measurements . Lcor  =  LEWIS_COR (patm) ; 

/*  I've  arbitrarily  cliosen  to  use  the  miscellanceous  category 
to  store  clothing  input  -  jrf  */ 

/*  All  lists  but  miscellaneous  deleted  for  v0.61  */ 
input_selections . clothing_type__idx  =0;  /*  miscellaneous  list  *! 

input_selections . clothing_idx  =0;  /*  MERCURY  INPUT  */ 

/*  select  the  table  */ 
datable  = 

clothing_list  [input__selections  .  clothing_type_idx]  .clo_lst; 

/*  load  the  table  w/  clothing  coefficients  */ 
datable  [  input___selections  .  clothing_idx]  .  itc  =  hsm__in.clo; 
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/*  clothing  velocity  coefficient  */ 
datable [input_selections . clothing_idx] . itvc  =  hsm_in.gc; 

/*  clothing  permeability  coefficient  */ 
datable [input_selections . clothing_idx] . imc  =  hsm_in.iclo; 

/*  clothing  permeability  velocity  coefficient  */ 
cl_table  [input__selections . clothing_idx]  . imvc  =  hsm_in.gi; 

/*  set  the  work  rate  index 

if  (hsm_in . workrate  <  150. OF)  resting 

input_selections . work_rate_idx  =  0; 

else  if  (hsm_in . workrate  <  250. OF)  v.  light 

input_selections . work_rate_idx  =  1 ; 
else  if  (hsm_in . workrate  <  425. OF)  light 

input_selections . work_rate_idx  =  2; 

else  if  (hsm_in. workrate  <  600. OF)  moderate 

input_selections . work_rate_idx  =  3; 

else 

input_selections . work_rate_idx  =  4;  heavy  */ 

if  (HSPredictionO  ==  FAILURE)  /*  indicate  FAILURE  to  MERCURY  */ 
max_work  =  -999; 

/*  write  hsdac.out  */ 
stream  =  fopen(fout,  "w"); 
if  (stream  ==  NULL)  { 

strcpy (buffer ,  fout) ; 
printfC'Can  not  open  %s\n" ,  buffer); 
exit (1)  ; 

} 


fprintf (stream,  "%4d\t\t!  Max  safe  work  time  [min] \n" ,  max_work) ; 
fprintf (stream,  "%4d\t\t!  Work-rest  cycle  [min] \n" ,  work_rest); 
fprintf (stream,  "%5.2f\t\t!  Water  rations  [canteens ] \n" ,  max_water) ; 
fprintf (stream,  "%5.1f\t\t!  Heat  strain  casualty  [%%]\n",  heat_cas) ; 
fprintf (stream,  "%5.2f\t\t!  Dewpoint  [C]\n",  measurements . dewpoint ) ; 
fprintf (stream,  "%5.2f\t\t!  Asp_wet_bulb  [C] \n" , measurements . awet_bulb) ; 
fprintf (stream,  "%5.2f\t\t!  Nat_wet_bulb  [C] \n" ,  measurements . wet_bulb)  ; 
fprintf (stream,  "%5.2f\t\t!  Globe  [C] \n" ,  measurements.globe); 
fprintf (stream,  "%5.2f\t\t!  WBGT  index  [C] \n" ,  measurements . wbgt ) ; 

*/ 

fprintf (stream,  "MxWk[min]\t  MxWtr  [qt/hr] \t  W/R  [min] \t  CyWtr 
[qt/hr] \t  RstWtr  [qt/hr] \t  RstWtrSh  [qt/hr] \tHeat  Cas  [%%] \t  Tdew 
[C]\t  Twba[C]\t  Twbn[C]\t  Tg  [C] \t  WBGT[C]\t  MWkRec  [min] \t 
MWkRecS [min] \t  patm\n") ; 

fprintf (stream,"%4d\t%5.1f\t%4d\t%5.1f\t%5.1f\t%5.1f\t%5.1f\t%5.1f\ 
t%5.1f\t%5.1f\  t%5.1f\t%5.1f\t%5.0f\t%5.0f\t%5.2f\n",  max__work, 
max_water,  work_rest,  water_req,  rest_water,  rest__water__sh, 
heat__cas ,  measurements,  dewpoint,  measurements . awet_bulb , 
measurements . wet_bulb,  measurements,  globe,  measurements . wbgt , 
maxwk__rcy,  maxwk__rcy__sh,  patm)  ; 

f  dose  ( stream)  ; 
exit ( 0 ) ; 
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#define  min(x,y)  ( (x  <  y)  ?  x  :  y)  /*  Return  the  minimum  of  2  values  */ 
#define  max(x,y)  ( (x  >  y)  ?  x  :  y)  /*  Return  the  maximum  of  2  values  */ 

static  float  xdata  Rload;  /*  Solar  factor  */ 

static  float  xdata  TmMWK_41;  /*  Max  work  for  light  casualties  */ 

/*********************************************************************** 
•*’*■*•  ★ 

*  Function  Name: 

*  HsPrediction 

* 

*  Description: 

*  This  is  the  main  HS  prediction  routine.  It  initializes  required 

*  variables  and  calls  the  routines  which  calculate  work/rest  and 

*  water  requirements. 

* 

*  Inputs : 

*  none 

* 

*  Outputs:  SUCCESS,  if  successful 

*  FAILURE,  if  not 

* 

*  *  *  *  / 

int  HSPrediction  () 

{ 

HSMALGO_STRUCT  xdata  hsm; 

ComputeAspiratedWetBulb {&measurements . awet_bulb) ; 

ComputeDewPoint () ; 

if  (ComputeNaturalWetBulb ( )  ==  FAILURE) 
return (  FAILURE  ); 

if  (ComputeGlobe {)  ==  FAILURE) 
return (  FAILURE  ) ; 

GetSkinTemp  (&hsm) ; 

CalcBaseValues  (ihsm) ; 

DoTheRest  (&hsm) ; 

return (  SUCCESS  )  ; 

}  /*  end  of  HsPrediction  routine  */ 

/i(**-ki(ici<*ici(iri(ic**i(*if*ici^*-kic**-k*ic-k-k-k-k-kic-kic-k-k-kic-k-k'k'k'k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k-k 
■k  -k  ic  -k  -k 

*  Function  Name: 
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DoTheRest 


■k 

*  Description: 

■k 

*  Inputs; 

*  none 

* 

*  Outputs: 

*  none 

* 


kkkkkkkkk-kkk-k-k-k'k-kk-kk'k'k-k-kk-k-k-k-kk-kk-k-k'k-k-k'k-kic-k-kif-k-k-k-k-k-k-kk-kk-kkkkkk-k-kic-kic-k-kk-k-k-k-k-k 


k  k  k  k  ! 

void  DoTheRest  (HSMALGO_STRUCT  *hsm) 

{ 


float 

float 

float 

float 

float 

float 


idata  temporary; 
idata  TmCWKold; 
idata  TmCWKnew; 
idata  TmMWKest; 
idata  TLest; 

MWTL ,  MCTL ; 


/*  localize  the  MAX_WORK_TEMP_LIMIT  and  MAX_CYCLIC_TEMP_LIMIT 
input  */ 

/*  these  were  hardwired  at  39.0  and  38.5  C  resp .  in  orginal  hsm_pc 

*/ 

MWTL  =  hsm_in . Tmxwork ; 

MCTL  =  hsm_in . Tmxcyc ; 

/****************************************************** 
*********************  *********** 
****■*••*■*■*•***★****■****★★★***•*••*•****★★*****★***'***★****•*•*★* 

/* 

*  Coefficients  for  Time  Series  Equations. 

*/ 

/*  time  delay  period  for  work  period.  */ 
hsm->TDWK  =  3480.  /  hsm_in . workrate ; 

/************************************************ 

*  only  for  initializing  the  first  work  period! 

************************************************y' 

hsm->TreoCWK  =  37.  +  (hsm->TrefRY  +  hsm->TrefARY  -  37.)  * 

(pow  (0.1,  pow  (0.4, ( (hsm->TDWK-30. ) /60 . ) ) ) ) ; 

/*  time  constant  for  the  work  period.  */ 

/*  vO . 6  change:  120  ->  225:  change  back  to  120,  May  2001-wtm  )  */ 

hsm->KWK  =  (1.  +  3.  *  exp  (0.3  *  (hsm- >TreoCWK  -  hsm- >Tref WK) ) )  / 

12  0  .  ; 


/*  rectal  temperature  difference  for  work  period.  */ 
hsm->DTreWK  =  hsm->TrefWK  +  hsm->TrefAWK  -  hsm- >TreoCWK; 

^kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk-kkkkkkkkkkkkkkkkkkkkkkkkk 
kkkkkkkkkkkkkkkkkkkkk  CgR  ******************** 
***************************************************************/ 
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/*  compute  time  delay  period  for  recovery  period  */ 
if  (hsm->CP  <  0 . ) 

hsm->TDRY  =  15.; 

else 

hsm->TDRY  =  15.  *  exp  (-0.5  *  hsm->CP); 

/*  time  constant  for  recovery  period  */ 

hsm->KRY  =  (1.  -  exp  (-1.5  *  fabs  (hsra->CP) ) ) /40 . ; 

/*************************************************************** 
*********************  mwk  ******************** 

***************************************************************y' 

TmMWK_41  =  MxWK  (MWTL,  hsm->DTreWK,  hsm- >TreoCWK,  hsm->TDWK, 

hsm->KWK) ; 

/*************************************************************** 
*********************  QY  ******************** 

***************************************************************y' 

/* 

*  Compute  the  work/rest  cycle. 

*/ 

hsm->TmCWK  =  0.0; 
hsm->TmCRY  =  0.0; 
if  (hsm->CP  >0.)  { 

TmCWKold  =  0 . ; 

TmCWKnew  =  0 .  ; 

if  (MCTL  >  hsm->TrefWK  +  hsm- >Tref AWK) 

/*  Initial  rectal  temp,  for  current  recovery  period  */ 
hsm->TreoCRY  =  hsm->TrefWK  +  hsm- >Tref AWK; 

else 

hsm->TreoCRY  =  MCTL; 

/*  rectal  temperature  difference  for  recovery  period  */ 
hsm->DTreRY  =  hsm->TrefRY  +  hsm->TrefARY  -  hsm- >TreoCRY; 

do  { 

/*  temporary  =  RECOVERY_TIME  in  MERCURY  ADA  -  jrf  */ 
temporary  =  hsm->TDRY  +  hsm->TDWK  +60.  -  0.5  * 

(TmCWKold  +  TmCWKnew) ; 

if  (hsm->TDRY  +5.  <  temporary) 

tmtry  =  temporary;  /*  !!  THIS  IS  TIME  -  jrf  */ 

else 

/*  initial  rectal  tmp  for  current  work  period  */ 
tmtry  =  hsm->TDRY  +  5.;  /*  !THIS  IS  TIME-  jrf  */ 

/*  NOTE:  hsm->TreoCWK  has  units  of  time  on  the  RHS  and 
temp  on  the  LHS  -  jrf  */ 
hsm->TreoCWK  =  TreRY  (tmtry,  hsm->TreoCRY, 
hsm->DTreRY,  hsm->KRY,  hsm->TDRY) ; 

hsm->DTreWK  =  hsm->TrefWK  +  hsm->TrefAWK  -  hsm- 
>TreoCWK; 
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TmMWKest  =  MxWK  (MCTL,  hsm->DTreWK,  hsm- >TreoCWK,  hsm- 
>TDWK,  hsm->KWK) ; 

TLest  =  TreWK( (TmMWKest  +  hsm->TDRY) ,  hsm- >TreoCWK, 
hsm->DTreWK,  hsm->KWK,  hsm->TDWK); 

TmCWKold  =  TmCWKnew; 

TmCWKnew  =  MxWK  ((MCTL,  -  0.5  *  (TLest  -  MCTL)), 

)ism- >DTreWK,  hsm- >TreoCWK, 
hsm->TDWK,  hsm->KWK) ; 

}  while  (  tabs  (TmCWKold-TmCWKnew)  >=  1.0  ) ; 

/*(  round  (TmCWKold,  0)  !=  round  (TmCWKnew,  0))*/ 

hsm->TmCWK  =  round  (TmCWKold,  0)  ; 
if  (10.  >60.  -  hsm->TmCWK) 

hsm->TmCRY  =  round  (hsm->TDRY  +5.,  0); 

else 

hsm->TmCRY  =60.  -  hsm->TmCWK; 

}  /*  end  if  CP  >  0  */ 

*  *  *  *  * 

•k  "k  ir  if  -k 

TmMWKest  =  TmMWK_41; 

if  (TmMWKest  >  300.0)  TmMWKest  =  300.0; 

if  (TmMWKest  >=  300.0)  { 

hsm->TmCWK  =  hsm->TmCRY  =  0  ; 
work_rest  =  60; 
max_work  =  300; 

}  else  if  (hsm->TmCWK  <  5)  { 

hsm->TmCWK  =  hsm->TmCRY  =  0; 
work_rest  =  0; 


if  (TmMWKest  > 

300 .0) 

max  work  = 

3  00; 

else 

max  work  = 

(int)  (round 

(TmMWKest, 

0)  )  ; 

}  else  { 

if  (hsm->TmCWK 

>=  300.0)  { 

hsm- >TmCWK 

=  hsm->TmCRY 

=  0; 

work  rest  = 

60; 

max  work  = 

300; 

}  else  { 

work  rest  = 

(int)  (round 

(hsm->TmCWK,  0) 

J 

i f  ( TmMWKe  s  t  > 

300  .  ) 

max  work  = 

300  ; 

else 

max  work  = 

} 

(int)  (round 

(TmMWKest , 

0)  )  ; 

water  req  =  max 

water  =  0.0; 

max_water  =  (round  {hsm->WtrWK,  1) ) ; 
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rest_water  =  (round  {hsm->WtrRY,  1) )  ; 

rest_water_sh  =  (round  (hsm->WtrRYSh,  1) ) ; 

if  (hsm->TmCWK  +  hsm->TmCRY  !=  0) 

water_req  =  round  ( (hsm- >TmCWK  *  hsm->WtrWK  + 
hsm->TmCRY  *  hsm->WtrRY)  /  (hsm->TmCWK  +  hsm->TmCRY)  ,  1)  ; 

}  /*  end  of  main  */ 

/*********************************************************************** 

*  Function  Marne: 

*  CalcBaseValues 

* 

*  Notes : 

*  Modified  4/13/98  jrf 

*  1)  added  solider  area  variable  (Asold)  and  passed  to  the  three 

*  functions  requiring  it:  Hrc,  Emax,  and  WTR. 

***********************************************************************y' 

void  CalcBaseValues  (HSMALGO_STRUCT  *hsm) 

{ 

#define  ALPHA  0 . 95F  /*  soldier  absorptivity  */ 

float  idata  Pa;  /*  Ambient  air  vapor  pressure  */ 

float  idata  HrcWK; 

float  idata  HrcRY; 

float  idata  EreqWK; 

float  idata  EreqRY; 

float  idata  EreqRYSh;  /*  NEW  wtm  Aug  01*/ 

float  idata  EmaxWK; 

float  idata  EmaxRY; 

float  idata  temp_wk; 

float  idata  temp_ry; 

float  idata  metab_w; 

float  Asold,  RloadSold, wt , Tmxwork, ereq ; 

Asold  =  DUBOIS (hsm_in . ht ,  hsm_in.wt); 
wt  =  hsm_in.wt; 

Tmxwork  =  hsm_in . Tmxwork ; 

RloadSold  =  Asold/1.8  *  Rload;  /*  Watts,  (Rload  is 

global  for  1.8  m2  man)  */ 

/*  measurements . Lcor  =  LEWIS_C0R (patm) ; */ 


/* 

*  Humidity  calculations. 

*  Pa=ambient  air  vapor  pressure  in  mmHg  (torr) 

*/ 

Pa  =  measurements . relhum*SVP (measurements . dry__bulb) / 10 0 . ; 

/* 

*  Get  the  metabolic  rate  in  watts  for  the  work  period. 
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*  (WTM  disabled  hard  categories : resting  =  105  very  light  =  150  light 
=  250  moderate  =  425  heavy  =  600) 

*/ 


metab_w  =  hsm_in . workrate ; 

/* 

*  Clothing  related  calculations. 

*/ 

temp_wk  =  Veff  (metab__w)  ; 
temp_ry  =  Veff(105.); 

/*  convective  and  radative  exchange  for  work  period.  */ 

/*  It {temp_wk)  is  clothing  insulation  value  for  work  period  */ 
HrcWK  =  Hrc(It(patm,  temp_wk) ,  hsm,  Asold) ; 

/*  convective  and  radative  exchange  for  recovery  period.  */ 

/*  It (temp_ry)  is  clothing  insulation  value  for  recovery  period  */ 
HrcRY  =  Hrc(It{patm,  temp_ry) ,  hsm,  Asold); 

/*  maximum  evaporative  loss  for  work  period.  */ 

/*  Cevap (temp_wk)  is  clothing  im/clo  value  for  work  period  */ 
EmaxWK  =  Emax(patm,  Cevap (temp_wk) ,  Pa,  hsm,  Asold); 

/*  maximum  evaporative  loss  for  recovery  period.  */ 

/*  Cevap (temp_ry)  is  clothing  im/clo  value  for  recovery  period  */ 
EmaxRY  =  Emax(patm,  Cevap { temp_ry) ,  Pa,  hsm,  Asold) ; 

temp_wk  =  U(temp_wk);  /*  clothing  efficiency  factor  for  work 

period  */ 

temp_ry  =  U(temp_ry);  /*  clothing  efficiency  factor  for 

recovery  period  */ 

/*  required  evaporative  loss  for  work  period.  */ 

EreqWK  =  Ereq (HrcWK,  metab_w,  temp_wk,  RloadSold) ; 

/*  required  evaporative  loss  for  recovery  period.  */ 

EreqRY  =  Ereq(HrcRY,  105.,  temp_ry,  RloadSold); 

/*  NEW:  required  evaporative  loss  for  rcovery  in  shade  */ 

EreqRYSh  =  Ereq (HrcRY,  105.,  temp_ry,  35.); 

/*  NEW:  recovery  time  from  max  work  in  the  sun  */ 

maxwk__rcy  =  MxWRKRY  (wt,Tmxwork,  EmaxRY,  EreqRY); 

/*  NEW:  recovery  time  from  max  work  in  the  shade*/ 

maxwk_rcy_sh  =  MxWRKRY  (wt,Tmxwork,  EmaxRY,  EreqRYSh); 
hsm~>CP  =  0.015  *  (EmaxRY  -  EreqRY) ;  /*  Cooling  power  */ 

hsm->CPSh  =  0.015  *  (EmaxRY  -  EreqRYSh) ;  /*  Cooling  power  in 
shade  */ 

/* 

*  Final  Rectal  and  Delta  Final  Rectal  for  Acclimatization. 

*  / 

/*  Final  rectal  temperature  for  work  period.  */ 
hsm->TrefWK  =  Tref (metab_w,  temp_wk,  HrcWK,  EreqWK,  EmaxWK, 
RloadSold) ; 
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/*  Final  rectal  temperature  for  recovery  period.  */ 
hsm->TrefRY  =  Tref(105.,  temp_ry,  HrcRY,  EreqRY,  EmaxRY, 

RloadSold)  ; 

/*  Final  rectal  temperature  add-on  for  work  period.  */ 
hsm->TrefAWK  =  DTref (hsm- >TrefWK,  EmaxWK) ; 

/*  Final  rectal  temperature  add-on  for  recovery  period.  */ 
hsm->TrefARY  =  DTref (hsm- >TrefRY,  EmaxRY); 

/* 

*  Water  Requirements . 

*/ 

/*  water  requirements  in  qt/hr  for  work  period.  */ 
hsm->WtrWK  =  WTR  (EreqWK,  EmaxWK,  Asold) ; 

/*  water  requirements  in  qt/hr  for  recovery  period.  */ 
hsm->WtrRY  =  WTR  (EreqRY,  EmaxRY,  Asold) ; 

/*  water  requirement  in  quarts  for  recovery  period  in  the  shade  */ 

hsm->WtrRYSh  =  WTR  (EreqRYSh, EmaxRY,  Asold); 

/*  compute  heat  casualty  */ 
heat_cas  =  heatcas (hsm- >TrefWK  +  hsm- >Tref AWK) ; 

}  /*  end  of  calc_base_values  */ 


/*********************************************************************** 

***** 

*  Function  Name: 

*  ComputeAspiratedWetBulb 

* 

*  Description:  Computes  wet  bulb  from  dry  bulb  and  RH. 

* 

*  Notes : 

*  1)  the  factor  2.0  is  the  Lewis  number.  It  has  units  of  degrees  C  / 
mm  Hg 

*  2)  Lewis  number  correction  for  atmospheric  pressure  Lcor=  Patm^(-.45) 

************************************************************************ 
*  *  *  *  y 

void  ComputeAspiratedWetBulb  (float  *Twb) 

{ 

float  idata  kl ; 
float  idata  k2 ; 
float  idata  delta; 
float  idata  awet__bulb; 

kl  =  measurements . Lcor  *  2.0  *  (SVP  (measurements . dry_bulb)  * 
measurements  .  relhum  /  100.0)  +  measurements  .  dry__bulb  ; 
k2  =  0.0; 
delta  =  10.0; 
awet_bulb  =  0 . ; 

while  (1)  { 

while  (kl  -  k2  >=  0 .  )  { 

k2  =  measurements .  Lcor  *  2.0  *  SVP  {awet_bulb)  +  awet___bulb; 
awet_bulb  +=  delta; 
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if  (awet_bulb  >  100.0)  break; 

} 

awet_bulb  -=  measurements . Lcor  *  2.0  *  delta; 

delta  /=  10.0; 

if  (delta  <  0.0001)  break; 

awet_bulb  -=  delta; 

k2  =  0.0; 

) 

measurements . awet_bulb  =  awet_bulb; 


void  ComputeDewPoint ( ) 

/******************************************************************** 
computes  the  dewpoint  from  the  dry  bulb  temperature  and  relative 
humidity 

***********************************************************************/ 

{ 

float  Pdp,  logPdp; 

Pdp  =  0.01  *  measurements . relhum  *  SVP (measurements . dry_bulb) ; 

/*  mm  Hg  */ 

logPdp  =  (float) logic (Pdp) ; 

measurements . dewpoint  =  (155. OF  -  235. OF  *  logPdp)  /  (logPdp  - 
8 . 1076F) ; 

} 


int  ComputeNaturalWetBulb ( ) 

/*********************************************************************** 
computes  the  natural  wet  bulb  temperature  using  the  model  of  Gonzalez 
et  al . 

Computes  by  solving  the  heat  balance 

f  =  Qin(Twb)  -  Qout (Twb)  =  0 
using  the  secant  iterative  method 

x_n+l  =  x__n  -  f  (x_n)  /  m 

m  =  (f  (x  n)  -  f  (x_n“l)  )  /  (x  n  -  x  n-1)  . 


where , 

Qin(Twb)  =  Qconv  +  Qrad 
Qout (Twb)  =  Qevap 

{ 

ttdefine  MAXITER2  10 
int  ict; 

float  Ta,  Trk,  dp,  hconv,  hevap,  Trk4 ,  Pdp; 
float  xnew,  xold,  fnew,  fold,  m,  reldif,  scale; 
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Ta  =  measurements . dry_bulb; 

Trk  =  measurements . mrt  +  273. 16F; 
dp  =  measurements . dewpoint; 

hconv  =  42.024F  *  pow { (patm*measurements .  wind_speed)  ,  0.466); 
hevap  =  measurements . Lcor  *  2.2F  *  hconv; 

Trk;4  =  Trk  *  Trk  *  Trk  *  Trk; 

Pdp  =  SVP (  dp  ) ;  / *  mm  Hg  * / 

/*********************************************************************** 

Let  our  first  guess  for  wet  bulb  be  equal  to  the  dewpoint. 

For  semantic  purposes  we  name  xnew  first  and  xold 
second  to  land  the  first  iteration  on  (xold,  fold) . 


***********************************************************************y' 

xnew  =  dp ; 

fnew  =  f_wetbulb (hconv,  hevap,  Ta,  Trk4 ,  Pdp,  xnew) ; 
xold  =  xnew  *  1.001; 

fold  =  f_wetbulb (hconv,  hevap,  Ta,  Trk4 ,  Pdp,  xold) ; 
m  =  (fold  -  fnew)  /  (xold  -  xnew) ; 

/*  need  to  define  "zero"  */ 

scale  =  fabs ( (xnew  +  273. 16F)  *  l.e-6); 

ict  =  0; 

lin_10 : 

if  (ict  >  MAXITER2) 

return (  FAILURE  ) ; 

xnew  =  xold  -  fold  /  m; 
reldif  =  fabs (xnew  -  xold) ; 
if  (reldif  >  scale)  { 

fnew  =  f_wetbulb (hconv,  hevap,  Ta,  Trk4 ,  Pdp,  xnew) ; 

m  =  (fnew  -  fold)  /  (xnew  -  xold) ; 

fold  =  fnew; 

xold  =  xnew; 

ict  =  ict  +  1; 

goto  lin_10; 

} 

measurements . wet_bulb  =  xnew; 
return)  SUCCESS  ) ; 

} 

/******************************************************************** 

f_wetbulb 

*******************************************************************y 

float  f_wetbulb (float  hconv,  float  hevap,  float  Ta,  float  Trk4 , 

float  Pdp,  float  Twb) 

/*********************************************************************** 
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the  heat  balance  equation  used  to  compute  the  natural  wet  bulb 
temperature 

WTM  changes  absorbtivity/emissicvity  of  wick  from  1  to  0.4*  STEFAN_B 

'*^*'*'************************************************ilf********-*'*****-)lr***-)ir 

{ 

float  Twbk,  Twbk4 ,  Pwb,  Qconv,  Qevap; 

Twbk  =  Twb  +  273. 16F; 

Twbk4  =  Twbk  *  Twbk  *  Twbk  *  Twbk; 

Pwb  =  SVP {  Twb  ) ; 

Qconv  =  hconv  *  (Ta  -  Twb)  +  0 . 4*STEFAN_B0LTZMAMN  *  {Trk4  - 
Twbk4)  ; 

Qevap  =  hevap  *  (Pwb  -  Pdp) ; 
return (Qconv  -  Qevap); 


/*********************************************************************** 

*  Function  Name: 

*  ComputeGlobe 

* 

*  Computes  the  globe  temperature  for  a  6  inch  globe  given  a  mean 

*  radiant  temperature  and  the  properties  of  the  air. 

* 

*  Also  computes  Rload  (a  global  variable),  the  radiant  load  [W/m2] 

*  on  the  solider,  and  the  wbgt  index. 

* 

*  Outputs:  SUCCESS,  if  successful 

*  FAILURE,  if  not 

* 

**********************************************************************^ 
int  ComputeGlobe  (void) 


float  d; 
globe  */ 
float  V; 

*/ 

float  xdata  Tr; 
temperature  */ 
float  xdata  Tgt ; 
float  xdata  Ta; 
temperature  */ 
float  xdata  Hcg; 
exchange  coef  */ 
float  xdata  Hrg; 
exchange  coef  */ 
float  Ta2 ,  Tr2 ; 

Ta  =  measurements . dry_bulb  +  2 
Ta2  =  Ta  *  Ta; 


/*  diameter  of 

/*  velocity  of  air 

/*  mean  radiant 

/*  globe  temperature  */ 
/*  ambient 

/*  convective  heat 

/*  radiative  heat 

3.16; 
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Tr  =  measurements . mrt  +  273.16; 

Tr2  =  Tr  *  Tr; 

d  =  0.1524;  /*  meters,  for  6"  globe  */ 

V  =  measurements . wind_speed;  /*  m/s  */ 

Hrg  =  0.95  *  STEFAN__BOLTZMANN  *  (Tr  +  Ta)  *  (Tr2  +  Ta2 )  ; 

Hcg  =  6.32  *  pow  (d,  -0.4)  *  pow  (  (patm  *  v) ,  0.5);  /* 

patm  is  bar  pres  correction  */ 

Tgt  =  (Hcg  *  Ta  +  Hrg  *  Tr) / (Hcg  +  Hrg)  ; 


measurements . globe  =  Tgt; 
measurements . globe  -=  273.16; 

/* 

*  Determine  the  radiant  load  [W]  using  USARIEM  Technical  Report 
No.TOl/13 
*/ 


Rload  =  -0.071  *  pow  ( (Tr  -  Ta) ,  2)  +  10.432  *  (Tr  -  Ta) ;  /*  W  */ 


/* 

*  Determine  the  WBGT. 
*/ 

measurements . wbgt  = 


0.7  *  measurements . wet_bulb  +  /*  outdoor  expr  */ 
0.2  *  measurements . globe  + 

0.1  *  measurements .dry_bulb; 


return (  SUCCESS  ); 

} 

/it********************************************************************** 
*  *  *  *  * 

*  Function  Name : 

*  SVP 

* 


*  *  ★  *  y 

float  SVP  (float  T) 

{ 

return  (pow  (10.,  (8 . 1076- (1750 . 286/ (T+235 .))))) ; 

}  /*  end  of  SVP  routine  */ 

/*********************************************************************** 
■k  -k  -k  -k  -k 

*  Function  Name: 

*  Veff 

* 

'kkk'k-k-k-k-kkkk-k-kkkkk-k-k-kkiik-k-k-k-k^-k-k-k-k-kifk-k-k-k-k-k-k-k-k'k-k-k-kk-k'kk-kkkkk'kkkkk-k-kk-kkk-kiikkk 
k  k  k  -k  ^ 

float  Veff  (float  M) 
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{ 

return  ( (measurements . wind_speed) +0 . 004*  (M-105 . ) )  ; 
}  /*  end  of  Veff  routine  */ 


/*********************************************************************** 

*  Function  Name: 

*  It 

* 

***********************************************************************/ 
float  It  (float  patm,  float  V) 

{ 

float  ItC; 
float  Itvc; 

clothing_table  *cl_table; 

/* 

*  Set  default  clothing  coefficients. 

*/ 

cl_table  = 

clothing_list  [input_selections . clothing_type_idx]  . clo_lst ; 

/*  clothing  coefficient  */ 

Itc  =  cl_table [input_selections . clothing_idx] . itc ; 

/*  clothing  velocity  coefficient  */ 

Itvc  =  cl_table [input_selections . clothing_idx] . itvc ; 

return  (Itc  *  pow{(patm  *  V),  Itvc) )  ; 

}  /*  end  of  It  routine  */ 

/*********************************************************************** 

*  Function  Name: 

*  Cevap 

* 

***********************************************************************/ 
float  Cevap  (float  V) 

{ 

float  Imc; 
float  Imvc; 

clothing_table  *cl_table; 

/* 

*  Set  default  clothing  coefficients. 

*/ 

datable  = 

clothing_list [input_selections . clothing_type_idx] .clo_lst; 

/*  clothing  permeability  coefficient  */ 

Imc  =  cl_table  [input__selections  .  clothing__idx]  .  imc  ; 

/*  clothing  permeability  velocity  coefficient  */ 

Imvc  =  cl__table  [input_selections  .  clothing___idx]  .  imvc  ; 
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return  {Imc*pow(V,  Imvc) ) ; 

}  /*  end  of  Cevap  routine  */ 

/*******************************************************************j.*** 

*  Function  Name: 

*  U 

* 

float  U  (float  V) 

{ 

float  ItC; 
float  Itvc; 

clothing_table  *cl_table; 

/* 

*  Set  default  clothing  coefficients. 

*/ 

cl_table  = 

clothing_list [input_selections . clothing_type_idx] . clo_lst ; 

/*  clothing  coefficient  */ 

Itc  =  cl_table [input_selections . clothing_idx] .itc; 

/*  clothing  velocity  coefficient  */ 

Itvc  =  cl_table [input_selections . clothing_idx] . itvc ; 

return  ( ( 0 . 41/Itc) *pow (V,  ( - ( 0 . 43+Itvc) ) ) ) ; 

}  /*  end  of  U  routine  */ 

/*********************************************************************** 

*  Function  Name: 

*  Hrc 

* 

*  Notes : 

*  Modified  4/13/98  jrf 

*  1)  the  original  constant  "11.772"  has  been  replaced  by  the  product 

*  6 . 4  5  *  A 

*  where 

*  A  is  the  DuBois  surface  area 

*  This  is  consistent  with  the  Ada  implementation  of  HSDAN . 

***********************************************************************^ 
float  Hrc  (float  It,  HSMALGO_STRUCT  *hsm,  float  A) 

{ 

return  (6.45F  *  A  * (measurements . dry_bulb  -  hsm- >Tskin) / It ) ; 

}  /*  end  of  Hrc  routine  */ 


36 


/*********************************************************************** 

*  Function  Name: 

*  Ereq 

* 

***********************************************************************/ 
float  Ereq  (float  Hrc ,  float  M,  float  U,  float  RloadSold) 

{ 

return  (Hrc  +  M  +  U  *  RloadSold) ; 

}  /*  end  of  Ereq  routine  */ 

/*********************************************************************** 

*  Function  Name: 

*  Emax 

* 

*  Notes : 

*  Modified  4/13/98  jrf 

*  1)  the  original  constant  "25.935"  has  been  replaced  by  the  product 

*  14 . 21  *  A 

*  where 

*  A  is  the  DuBois  surface  area 

*  This  is  consistent  with  the  Ada  implementation  of  HSDAN. 

*  2)  the  constant  45.8  has  been  replaced  by  the  function  SVP,  given  the 

*  skin  temperature. 

*  3)  Emax  factor  for  high  terrestrial  altitude,  patm  ^-0.45  wtm  Nov 
2001 

***********************************************************************y 
float  Emax  (float  patm,  float  Cevap,  float  Pa,  HSMALGO_STRUCT  *hsm, 
float  A) 

{ 

float  result; 


result  =  (pow(patm,  -0.45)  *  14.21F  *  A  *  Cevap  *  (SVP  (hsm->Tskin) 

-  Pa) )  ; 

/*  result  =  (14.21F  *  A  *  Cevap  *  (SVP  (hsm->Tskin)  -  Pa));*/ 
return (  result  ); 

/*  4/13/98  jrf  return  (25.935  *  Cevap  *  (45.8  ~  Pa));  */ 

/*  3/19/97  KSH  return  (25.935  *  Cevap  *  (SVP  (hsm->Tskin)  -  Pa));  */ 

}  /*  end  of  Emax  routine  */ 

/********************************************************************** 

*  Function  Name: 

*  Tref 

* 

***********************************************************************^ 
float  Tref  (float  M,  float  U,  float  Hrc,  float  Ereq,  float  Emax, 

float  RloadSold) 
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return  (36.75  +  0.004*M  +  0.0025*  U  *  RloadSold  +  0.0011*Hrc  + 
0 . 8*exp ( 0 . 0047* (Ereq  -  Emax) ) ) ; 

}  /*  end  of  Tref  routine  */ 

/******************************************************************* 

*  Function  Name: 

*  DTref 

* 

*  Notes : 

*  Modified  4/13/98  by  jrf 

*  1)  added  logic  for  Kdehyd 

*  2)  acclimation  assigned  directly  to  days  in  heat  (DIH) 

*  3)  broke  big  eq  into  smaller,  more  readable  pieces 

*  4)  removed 

*  if  (DTref_tmpry  <  O.OOOIF) 

*  DTref_tmpry  =  O.OF; 

*  logic 


float  DTref  (float  Tref,  float  Emax) 


float  DTref_tmpry; 

float  DIH; 

float  Kdehyd,  alp; 

/*  days  in  heat  */ 
DIH  =  hsm_in . acclim; 


/*  dehydration  factor  */ 
if  (hsm_in . dehyd  >=  5 . OF) 

Kdehyd  =  l.OF; 
else 

Kdehyd  =  0.2F  *  hsm_in . dehyd; 

alp  =  -0.3F  *  DIH  *  (1  -  Kdehyd); 

DTref_tmpry  =  0.0; 
if  (Emax  >  0 . )  ( 

DTref _tmpry  =  (0.5F  + 

1 . 2F  * 

(l.OF  -  exp  (0 . 5*  (37 . 15  -  Tref)))  * 
(l.OF  -  exp ( - 0 .005*  Emax ) )  )  ; 

DTref__tmpry  *=  (float)  exp  (alp)  ; 

DTref_tmpry  +=  0.1735F  *  hsm_in . dehyd  -  0.215F; 


/* 


if  (DTref_tmpry  <  O.OOOIF) 

DTref_tmpry  =  O.OF;  */ 


return  (DTref _tmpry) ; 

}  /*  end  of  DTref  routine  */ 


/***************************************.*******************J.******JJ****JJ 

*  Function  Name : 

*  WTR 

* 

*  Notes : 

*  Modified  4/13/98  jrf 

*  1)  the  original  constant  "50.921"  has  been  replaced  by  the  product 

*  27 . 9  *  A 

*  Modified  March  2001  wtm 

*  1)  original  Shapiro  Ereq  and  Emax  are  per  unit  BSA! 

*  where 

*  A  is  the  DuBois  surface  area 

*  This  is  consistent  with  the  Ada  implementation  of  HSDAN. 

***************************************************************.********y 

float  WTR  (float  Ereq,  float  Emax,  float  A) 

{ 

float  sweat; 

/*  new  variables,  body  surface  area-normalized  values  of  Ereq  and 
Emaxfor  Shapiro  Eq.  wtm  */ 

float  Ereqs ; 

float  Emaxs; 

Ereqs=Ereq/A; 

Emaxs  =  Emax/ A ; 

if  (Emaxs  <=  0.0) 

sweat  =  2000  .  ; 

else 

sweat  =  27.9F  *  A  *  Ereqs  *  pow (Emaxs,  -0.455); 

/* 

*  Limit  sweat  to  the  range  150  -  2000. 

*/ 

sweat  =  max  ((min  (sweat,  2000.)),  150.); 

/* 

*  Return  the  water  requirements . 

*/ 

return  (sweat  *  .0010567);  /*  liters  ->  canteens  (i.e.  qts)  */ 

}  /*  end  of  WTR  routine  */ 

/*********************************************************************** 

*  Function  Name: 

*  MxWK 

* 

***********************************************************************/ 
float  MxWK  (float  TL,  float  DTreWK,  float  TreoCWK,  float  TDWK,  float 
KWK) 
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float  temp; 


if  (DTreWK  >0.)  { 

temp  =  (DTreWK  +  TreoCWK  -  TL)  /  DTreWK; 
if  (temp  >  0 . ) 

temp  =  log  (temp) ; 

else 

temp  =  -227.9;  /*  same  as  log  (16-99)  */ 

}  else  { 

temp  =  -227.9;  /*  same  as  log  (le-99)  */ 

} 

return  (TDWK  -  temp/KWK) ; 

}  /*  end  of  MxWK  routine  */ 


*  Function  Name: 

*  round 

•k 

kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk-k-k-k-k-kk-k-k-k-k-kk-k-k-kk-k-kk-k-k-kk-k-k-kk-kkk-k-k-k-kk-k-k/ 

float  round  (float  x,  unsigned  char  digit) 

{ 

X  *=  pow  (10.0,  (float) digit) ; 

X  =  (float) ({int)(x  +  0.5)); 

X  /=  pow  (10.0,  (float) digit) ; 

return  (x) ; 

}  /*  end  of  round  routine  */ 

/*********************************************************************** 

*  Function  Name: 

*  TreWK 

* 

***********************************************************************y' 

float  TreWK  (float  t,  float  TreoCWK,  float  DTreWK,  float  KWK,  float 
TDWK) 

{ 

return  (TreoCWK  +  DTreWK  *  (1.0  -  exp  (-KWK  *  (t  -  TDWK) ) ) ) ; 

}  /*  end  of  TreWK  routine  */ 

*  Function  Name: 

*  TreRY 

* 

kkkkkkkkkkkkkkkkkkkkkkkk-k-kk-k-k-k-k-k-k-kk-k-k-kk-kickkk-kk-kkkkkk'kk-k'kk-k-k-k-k-k-k-k-k-k'k-k-k-k-k/ 
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float  TreRY  (float  t,  float  TreoCRY,  float  DTreRY,  float  KRY,  float 
TDRY) 

{ 

return  (TreoCRY  +  DTreRY  *  (1.0  -  exp  {-KRY  *  (t  -  TDRY) ) ) ) ; 

}  /*  end  of  TreRY  routine  */ 


/****************************************************** 

*  Function  Name : 

*  MxWKRY 

* 


■k-k-k-k-k-k-k-k-k-k-k-k-k-k'k-k-k 


*  note: This  is  Lee's  StrTm, cooling  based  on  heat  content, 

*  0.965  is  0.83  body  sp  heat  /  0.8606  Kcal/hr/Watt 

*  wtm  Aug  2001 

* 


**********************-k-i:*-k-t:-k-k-k**-k*-lt***ic***-tc**-k****-ki,*-k****i,.)!.i,i,i,i,i,i^^^,^,i,^ 

float  MxWRKRY  (float  wt,  float  Tmxwork,  float  EmaxRY,  float  ereq) 


{ 


float  temp; 

temp  =  (0.95  *  wt  *  (Tmxwork  -  37 . 0 )/ (EmaxRY  -  ereq) ) *60; 
if  (temp  <  0  )  temp  =  0; 
return  (round (temp,  1) ) ; 


}  /*  end  of  MxWRKRY 


/*******************ir****-k-t!*-k-k********-t!**-H:-)!****-k*ir**ir-lri,i,i^i,i,.t:^.lr^.^.^.i.^^^^^ 

*  Function  Name: 

*  GetSkinTemp 


*  note:  dead  code 

*  Hold  Tskin  to  36.5 

void  GetSkinTemp  (HSMALGO  STRUCT  *hsm) 


#if  0 

float  TO; 
float  Hc; 


He  =  8.6  *  pow  (measurements . wind_speed,  0.5); 

TO  =  (4.7  *  MRT  +  Hc  *  measurements . dry_bulb)  /  (4.7  +  Hc) ; 

hsm->Tskin  =  25.8  +  0.267  *  TO ; 

/* 

*  Constrain  to  34  -  37. 

*/ 

hsm->Tskin  =  min  (37.0,  max  (34.0,  hsm- >Tskin) ) ; 
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#else 


hsm~>Tsk;in  =  3  6.5; 

#endif 

)  /*  end  of  GetSkinTemp  routine  */ 
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